Asymmetric Dirac cones in monatomic hexagonal lattices. 
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The study of nanostructures has contributed to the advance of an interdisciphnary science as 
nanotechnology is. Among those ones, graphene has been distinguished in the last years by its 
interesting properties and specially, in our own interest, the presence of Dirac cones in electronic 
dispersion relation. This material is considered the cornerstone in further scientific and techno- 
logical innovation. In this theoretical paper, crystallographic phase systems similar to graphene, 
such as silicon hexagonal monolayer (h-Si) or germanium (h-Ge), are analyzed, using the density 
functional theory (DFT), implemented in the code SIESTA, in the generalized gradient approxi- 
mation (GGA), into the Perdew-Burke-Ernzerhof functional (PBE) and optimized norm-conserving 
pseudopotentials. The results found permit us to report the chemical stability of h-Ge. Also, the 
lattice parameters, electronic dispersion relations and the density of states (DOS) for graphene, h-Si 
and h-Ge are reported. The existence of Dirac cones is seen in the electronic dispersion relation for 
each one of the studied hexagonal monolayers. The cones show lateral asymmetry in function of the 
direction around the K point. In particular, Fermi velocities are calculated for holes and electrons 
in function of the direction. 
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I. INTRODUCTION 

Currently, the nanotechnology boom as multidisci- 
plinary science is owing to the construction and charac- 
terization of a wide variety of materials with low dimen- 
sionality as much as organic, inorganic or mixed types^. 
Recently one of the most important results has been 
isolating the graphene through mechanical exfoliatiorP. 
This monolayer of carbon atoms prepared in hexago- 
nal form accounts for two dimensional structure material 
which existence should be impossible due to thermody- 
namic instability showed by L. Landau and R.Peierls' 
theor}0El. 

Graphene has been investigated either in theory or ex- 
perimentally wa}'^^^, finding in its dispersion relation, 
the absence of an energy gap, with an innovative char- 
acteristic such as the linearity of conduction and valence 
bands around Fermi level that is found in the reciprocal 
lattice i^-point. In these zones known as Dirac cones, 
electrons and holes behave as cuasiparticles without mass 
and with an effective speed c* ^ 10^ m/s, independent 
of its wavevector. In front of this important discovering, 
we wonder ourselves if other stable materials can exist in 
two dimensional form and show linearity in its electronic 
dispersion relation. 

Searching for those answers has motivated the devel- 
opment of this work to explore structural and electronic 
properties of materials which through years have had 
a huge diversity on technological applications such as: 
Si (Silicon) and Ge (Germanium), but now in hexag- 
onal monolayer forms (h-Si, h-Ge). In theoretical in- 
vestigation of low dimensional systems, formalisms of 
first principles or ah initio have been used in its Den- 
sity Functional Theory (DFT) formulation. In particular, 
the SIESTA code implementation {Spanish Initiative for 
Electronic Simulations with Thousands of Atoms ]^ is em- 
ployed in the academic environment^ as in the industriaP 



or applied investigatioiP^. This article is organized as fol- 
low. In section [n| some computer details relate with the 



method are introduced, in section |III| some relevant phys- 
ical parameters are reported. In sections |IV| and |V| we 
show the electronic dispersion relations. Density of States 
(DOS) and the analysis of the Dirac cones for these ma- 
terials. 



II. COMPUTER DETAILS 

In the study of electronic and structural properties of 
h-C, h-Si and h-Ge monolayers, we use SIESTA code, 
which useful to realize calculations of electronic struc- 
ture or vibrational properties and geometrical optimiza- 
tion by simulations of molecular dynamic. Besides, It 
employs efficient algorithms permitting a medium level 
computer cost, due to the time of linearly scale calcu- 
lus with the number of system atoms. In the study 
of our systems, the only initial information provided 
to the code is the number and the type of atoms that 
compounds each system with their approximated posi- 
tions in the hexagonal structure and lattice parameter. 
In the calculus of interchange and correlation energy 
is used the gradient generalized approximation (GGA) 
with the Perdew-Burke-Ernzerhof functional (PBE)P2. 
Valence electrons are described for localized numerical 
atomic orbitals (NAO) with a double— (" base polariza- 
tion (DZP) and core electrons are treated implicitly with 
norm- conserving pseudopotentials^^'^ and optimized cut 
off radii Tc^^ ^^. The analysis of the chemical stability of 
the structures is realized through a relaxation process, 
in which the unit cell atoms are moved slowly, step by 
step, minimizing the strength among them and with 200 
interactions per step. A chemical structure is considered 
stable for a maximum atomic force tolerance of 40 meV 
A^. In the analyzed systems, stability was found before 



2000 steps. 



III. LATTICE PARAMETERS 

Lattice parameters for different bidimensional hexag- 
onal sheets are found through a relaxation process by 
dynamic molecular, in which chemical stability is deter- 
mined of these systems according to values mentioned 
above. In particular, the data found for bond length ao, 
lattice constant a and maximum atomic force tolerance 
ftoi in the process of relaxation are reported in the ta- 
ble |Tj As a way of method validation are included the 
values found for graphene, which are in good agreement 
with previous theoretical and experimental reports up to 
999^31111. 



TABLE L Bond length ao, lattice constant a and maximum 
atomic force tolerance ftoi of hexagonal monolayers, including 
spin effects (cs) and absence of these (ss). 



is a = 3.3 A^i and a = 3.65 A for (0001) - ZrB'P^, 
When these results are compared with ours, we found 
differences of 0.46%, 14.55% and 5.49%, respectively. 



IV. DISPERSION RELATIONS AND DENSITY 
OF STATES 

Silicon (Si) and Germanium (Ge) structures in their 
three-dimensional form are equivalent to the diamond 
one with indirect energy gaps Eg of 1.12 eV and 0.67 eV 
to T = 300 ^33,24^ respectively, showing a semiconduc- 
tor behavior. For bidimensional materials, i. e. hexag- 
onal sheets, after the relaxation process which permits 
to obtain the lattice parameter, the electronic dispersion 
relation is calculated, either with spin effects (ss) or in- 
cluding them (cs). 

In particular, for graphene, spin effects are evidents to 
high energies while for values around Fermi levels, signif- 
icant changes are not observed in the DOS, see Figs. [l]y 
[2] Nevertheless, including spin effects generates relevant 
changes in the dispersion relation and specially in the 
Dirac cones, such as it can be observed in Fig. [3) where 
there are an enlargement for electronic dispersion rela- 
tions around K-point without including spin effects (ss), 
left panel, and including them (cs), right panel. Main 
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FIG. 1. Electronic dispersion relation and density of states 
(DOS) for a hexagonal monolayer of graphene without spin 
effects (ss). 



Monocapa ao (A) a (A) 


ftoi (eVA-1) 




h-C (ss) 1.419 2.458 


0.9x10"^ 




h-C (cs) 1.419 2.458 


0.002 




h-Si (ss) 2.230 3.860 


0.3x10"^ 


> 


h-Si (cs) 2.230 3.860 


0.2x10"^ 


w 


h-Ge (ss) 2.310 4.001 


0.001 


0£ 


h-Ge (cs) 2.310 4.001 


0.030 




In case of h-Si, the values reported depend on substrate 


w 


used in their growth. In substrates of A 


0(110) the lattice 




parameter reported is 3.88 A^^, using Ag {111) the value 
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FIG. 2. Electronic dispersion relation and density of states 
(DOS) for a hexagonal monolayer of graphene with spin effects 
(cs). 



observed changes refer to the Fermi speed both for elec- 
trons and holes and an lateral asymmetry of each one of 
the Dirac cones as function of approaching direction to 
K point. Previous theoretical results using models type 
tight binding^, do not show this lateral asymmetry, nei- 
ther the difference among speed values of the electrons 
and holes, which has been observed experimentally such 
as Kai-Chieh Chuang et a?^ and Jiam Xue et a?^ report. 
In the case of h-Si and h-Ge, a semimetal character- 
istic is observed, i. e., the non existence of an energy 
gap between valence and conduction bands, see Figs. [4] 
and |5j These ones get intersection punctually around 
K^ equivalent to what occurs in graphene. Characteris- 
tic that is confirmed in DOS for each case, as it can be 
seen in the right panel of the corresponding dispersion 
relations. Previous results report that h-Si is a direct 
bandgap semiconductor as show the table [llj Our results 
indicate that h-Si and h-Ge are semimetal systems, in- 
stead of spin effects are included or not. 









TABLE II. Energy gap values in eV. 






System 


Ref- 


Ref.29 


Ref.2° Ref.31 Ref.^^ Ref.^^ 


Ref.^^ 


Here 


h-Si 
h-Ge 


Semimetal 
Semimetal 


Semimetal 
Metal 


Metal Metal 0.064(Direct) 
Metal - Eg^O -0.444(Semimetal) 


Semimetal 
Semimetal 


Semimetal 
Semimetal 




W-0,5 



FIG. 3. Enlargement for electronic dispersion relation around 
i^-point for graphene, without including spin effects (ss), left 
panel, and including them (cs), right panel. 



Dispersion relation of h-Si and h-Ge show a hnear be- 
havior around Fermi level, giving rise to Dirac cones, 
FigsJ6] y [7| equivalent to results for graphene i. e., in 
these materials quantum electrodynamics experiments 
(QED) could be developed, also. 



> 
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FIG. 4. Electronic dispersion relation and density of states 
(DOS) for h-Si with spin effects (cs). 

In Figs. [6] and [7] there are enlargements of electronic 
dispersion relations around if-point for h-Si and h-Ge 
where is possible to visualize lateral asymmetry for the 
Dirac cones and different velocities between electrons and 
holes, without considering spin, left panel, and including 
its effect, right panel. These lateral asymmetries and 
differences are similar to the ones discussed for graphene. 

A quantitative analysis of dispersion relation at Dirac 
cones permits to determinate Fermi speeds for electrons 
and holes, in each of the studied materials. For this 
purpose, a linear regression is done in T ^ K and of 
M ^ K directions. Results are showed in the table \lT\\ 
where lateral asymmetry for Fermi speed is evidenced 
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FIG. 5. Electronic dispersion relation and density of states 
(DOS) for a Ge hexagonal monolayer with spin effects (cs). 



either for electrons or holes. The relative variation of 
Fermi speed to the electrons, in function of the wavevec- 
tor k = kxi -\- kyj^ between the two analyzed directions 
is 61,4% for h-Si(cs), of 45.9% for h-Ge(cs) and of 36.6% 
for graphene (cs). Equivalent differences are found for 
Dirac cones holes. In h-Si(cs) there is a relative variation 
of 57.5%, in h-Ge(cs) of 48.5% and of 28.7% for graphene 
(cs). 
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FIG. 6. Enlargement for electronic dispersion relation around 
i^-point for h-Si, without including spin effects (ss), left panel, 
and including them (cs), right panel. 

Experimental results report Fermi speeds of electrons 
in graphene of 1.093 x 10^ m/s^^, 0.79 x 10^ m/J^, 1.10 x 
1q6 ^/jnu^ 1,16 ± 0.01 X 10^ m/s and for holes of 0.94 ± 
0.02 X 10^ m/ J^. In h-Si experimental results give a value 
of 1.3 X 10^ m/s for Fermi speed of electrons and holes^. 
For h-Ge, so far we know, there are no experimental speed 
reports. 

The relative variation of Fermi speed for electrons and 
holes, in function of the /c-space directions, indicate that 



TABLE III. Fermi speeds Vf (xlO^m/s) for electrons and holes, in graphene, h-Si and h-Ge. Without spin effects (ss) and 
including them (cs). 



electrons 


holes 




Vfir^K) 


Vf{K^ M) 


Vfir^K) 


Vf{K^ M) 


ss 
h-C 

cs 


1.406 
1.365 


0.998 
0.999 


1.409 
1.385 


0.996 
1.076 


ss 
h-Si 

cs 


1.478 
1.509 


1.078 
0.935 


1.550 
1.533 


1.062 
0.973 


ss 
h-Ge 

cs 


1.644 
1.746 


1.047 
1.197 


1.569 
1.639 


0.981 
1.104 



this is not constant, presenting a spacial dependence or 
anisotropy. In recent experimental reports^^ 38 40 jg ]^gg^_ 
sured a renormalization of Fermi speed. This feature 
is not unique of the graphene, being shared by other 
monoatomic bidimensional hexagonal lattices such as: 
h-Si and el h-Ge. Theoretical results, not show here, 
indicate that diatomic hexagonal monolayers have not 
present this behavioi^^. 
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FIG. 8. Dirac cones in 3D for holes and electrons in h-Si and 
h-Ge with spin effects (cs). Color online 



FIG. 7. Enlargement for electronic dispersion relation around 
i^-point for h-Ge, without including spin effects (ss), left 
panel, and including them (cs), right panel. 



V. CONCLUSIONS 

Lateral asymmetry or spatial dependence of Fermi 
speed for electrons and holes is reflected in asymmet- 
ric Dirac cones, as it can been observed in Figs. IsFa) 
andpfb). These cones are obtained from DFT calculus 
using SIESTA implementation. In particular, we use the 
dispersion relations from some M, F and iC-points neigh- 
bors to i^-point at (|, ^y^). Fig 8[a) corresponds to h-Si 
case while Figjsjb) for h-Ge. In both cases, it is possible 
to identify two type of asymmetries, i) A lateral asymme- 
try that deform the Dirac cone far away from a circular 
shaped, ii) An an asymmetry electron- hole i. e., different 
Dirac velocities for valence and conduction bands. Hole 
Dirac cones show less height compared with the electrons 
ones, being this difference stronger for h-Ge, indicating 
that behavior type no massive Dirac Fermions for holes is 
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0X1 
0.2 



FIG. 9. Contours levels for the energy of for h-Si and h-Ge 
with spin effects (cs). Color online 



lost rapidly at the same time that the wavevector get dis- 
tant from K point. This characteristic can be important 



in experiments that involve magnetic fields or phonon as- 
sisted transitions. In Figs. [9] there are energy level curves 
or equi-energy contours at Dirac cones for holes, upper 
panels, and for electrons, lower panels, in h-Si Fig|9^ and 
h-Ge Fig. ^p. The intensity of colors represents the en- 
ergy magnitude in relation with the Fermi level. In each 
case 10 contour lines are shown. Clearly, there are not 
circular profiles respect to point K. Lateral asymmetry 
is represented for the wavevector magnitude variations, 
or the distance from K -point to each one level curves. 
In particular, Dirac cones radial asymmetry for holes is 
bigger in comparison with Electron Dirac cones. From 



Dirac equation for no massive Fermions given by^ 

E±{k) ^ ±VF\\k\\, 



na 



(1) 



wavevector variations over each one of these constant en- 
ergy lines implies Fermi speed must do it too, revealing 
a spatial asymmetry of this speed. 
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